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^ '. Abstract 

The local work function of a surface determines the spatial decay of the charge density at the 
Fermi level normal to the surface. Here, we present a method that enables simultaneous mea- 
surements of local work function and tip-sample forces. A combined dynamic scanning tunneling 
microscope and atomic force microscope is used to measure the tunneling current between an os- 
dilating tip and the sample in real time as a function of the cantilever's deflection. Atomically 
resolved work function measurements on a silicon (lll)-(7 x 7) surface are presented and related 
to concurrently recorded tunneling current- and force- measurements. 
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When two metallic electrodes are separated by a small vacuum gap and a bias voltage V t 



is applied, a quantum mechanical tunneling current flows l|, SQl- For metallic electrodes, 
the tunneling current increases approximately by a factor of 10 for each distance reduction 
of 100 pm. This sharp distance dependence is key to the atomic resolution capability of the 
scanning tunneling microscope (STM) |4|. In one dimension, the decay of the tunneling 
current with distance is roughly given by 

I = J exp (-2kz) (1) 

with 

k = ^2m$/h (2) 

where m is the mass of an electron and $ is the local work function or barrier height 2] . As 
Binnig has shown early in the development of STM, this exponential dependence holds over 
at least 0.5 nm or so for tunneling currents ranging from a few hundred pA to a few fiA 
Due to the atomic structure of matter, k is also a function of the lateral positions x and y 
as well as z (see Fig. 1) for very small tunneling distances. 

A straightforward method for measuring k(x, y, z) is to stop the lateral scan over a specific 
atom position and perform a I(z) measurement. However, this method is time consuming 
and prone to errors from creep and drift of the piezoelectric scanner. Pethica et al. 5J] 
have extended 'current-imaging-tunneling-spectroscopy'(CITS) by Hamers et al. an ac- 
method for current vs. voltage spectroscopy where a small ac-voltage is added to Vt to 
recover the density of states. In Pethica et al.'s method, the tip is oscillated at a fixed 
frequency / (typically on the order of a few kHz) according to x = x~q + Acos(2rcft) [5J. For 
this situation, Eq. [T] needs to be generalized to 

I = Jq exp (— 2kx). (3) 

The natural logarithm of the normalized current is then given by 

ln(I/I ) = -2kx - 2/L4cos(27r/t). (4) 

In general, we define an apparent decay 'constant'/? = r\j (2A), where r\ is the AC-component 
of ln(/// )- For the case described by Eq. |H 77 = 2kA. 



Atomic force microscopy (AFM) has progressed rapidly in the past years 
combination with other techniques like STM |9j or Kelvin Probe Mircoscopy 
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FIG. 1: Schematic view of In (I/Iq) above two silicon adatoms located at z = and z = 0.77 nm of a 
silicon (lll)-(7x 7) surface. The surface is oriented perpendicular to the z-axis. The current distri- 
)ution between each of the silicon states and a s-tip is approximately given by I oc (z/r) 2 exp (— 2nr) 
3]. The amplitude of the vertical oscillation is A. <p is the angle between the oscillation direction 
and k. For decreasing cos (<p), the decay strength of the current in vertical direction also decreases 
(see text). Note that generally k is not parallel to Vln (I/Iq) because of the angular factor of the 
wave function. The distance of the contour lines corresponds to a current increase of 1.82. 



feasible. Here, we combine AFM with Pethica's method by using a qPlus sensor [11| where 
the STM tip is mounted on a vibrating cantilever. The oscillation frequency is no longer 
fixed in this case, but varies by a frequency shift Af as determined by the atomic forces 
acting between tip and sample. The frequency shift data can then be related to forces 0, 
and simultaneous measurements of forces and decay constants are possible. 

Because the atomically resolved STM image is influenced by the atomic and subatomic 



structure of the tip and sample wavefunctions 



DO 



12 1, these structures also influence 



the decay constant images |5|. The tunneling current can be calculated with a plane- wave 



expansion of the sur: 
atomic orbitals 



121. 



ace wavefunctions ultimately depending on the superposition of 



For the atomic basis functions at the Fermi level, the exponential radial factor is 
exp (— Kr), where k = nrj r and r is the position vector with respect to the nucleus. When the 
ip moves by a, the corresponding conductivity distribution changes according to exp (— 2ko) 
2j]. Fig. ^ shows a cross-section of a model distribution of ln(/// ) for two atomic silicon 
states of a silicon (lll)-(7 x 7) surface, with the nuclei on the x-axis at z = 0, separated by 
0.77nm. The current dependence / oc (z/r) 2 exp (— Ikt) is a good fit for tunneling between 
a tip (s-tip) and each of the sp 3 -adatom states 2|. The tip is assumed to oscillate verti- 
cally, ip is the angle between the oscillation direction and k. For cos(v?) — > 1, r\ approaches 
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I/V logarithmic 
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FIG. 2: Principle of dynamic workfunction measurement. The tip is mounted on a quartz cantilever 
(not shown here, see [h| for details) with a stiffness of 1800 N/m which oscillates at a fixed 
amplitude A ~ 0.1 nm. The unperturbed resonance frequency is /o ~ 20 kHz. Due to the sinusoidal 
oscillation of the cantilever, the tunneling current is strongly modulated at a frequency / ~ /o (see 
Fig. 3). 

2KA = IkA in this model. Therefore, an apparent barrier height $ app = h 2 k 2 /(2m) can be 
extracted from the value of k at the position of the single atom. For decreasing cos(<y?), k also 
decreases in the theoretical model, since A is no longer parallel to the direction of the fastest 
decay of the current. In the model, odd higher powers of [ln(/// )](g) = [ln(J / Io)](xo+qA/ A) 
influence the value of well. 

In this publication, we present dynamic STM/AFM measurements performed in ul- 
trahigh vacuum at a pressure of p ~ 10~ 8 Pa and ambient temperature T ps 300 K. A silicon 
lll)-(7 x 7) surface is probed by a tungsten tip that is mounted on a qPlus force sensor 



11 1 . The qPlus sensor used here has a stiffness of k = 1800 N/m and an eigenfrequency 



of 15487 Hz. The oscillation amplitude of the cantilever needs to be approximately below 
0.1 nm such that / is above the noise level of the current amplifier even when the tip is at 
the far point of the surface. Operation at such small amplitudes poses a challenge as ampli- 
tude fluctuations tend to increase with decreasing amplitude, and only the use of cantilevers 
with a very large stiffness decreases amplitude fluctuations to an acceptable level ■8]- The 
experimental setup is shown in Fig. El 

Fig. |3] shows a simulation of the excursion of the tip mounted on a qPlus cantilever, 
the normalized tunneling current I and the natural logarithm of I for an ideal exponential 




FIG. 3: Tip deflection divided by amplitude, current divided by maximum current and natural 
logarithm of I/Io as a function of time t for t = to t = 1//. If the oscillation amplitude is large 
compared to the decay length of the tunneling current, the current appears as a series of Gaussian 
peaks. 

current- distance dependence. 

Fig. E] shows a simultaneous measurement of topography (constant average current), k, 
frequency shift Af and damping signal AE C per oscillation cycle. The angle between A and 
the surface normal was 9 ~ 20°. The cantilever-amplitude was A = 100pm±4pm. At a 
single adatom defect site in the right bottom of Fig. 0] (c), the frequency shi 
because of a very small tip-sample distance, where repulsive forces are acting (9, 



t is positive 

IMM- The 



presence of a single atomic defect proves that the tip has a single front atom and multitip 



effects can be ruled out. The average damping signal is AE C =170 meV and exceeds the 
internal damping of the cantilever by ~ lOmeV. On an atomic scale, k shows a strong 
variation. The maximum value of k is ~ 65% of the nominal value of k for this tip-sample 
system (workfunction $ fa 4.5 eV, [14]). The average of k is less than 50% of the nominal 
value of k. This lowering of the measured average R can not be explained with the hypothetic 
geometry factor cos(6 l ) for the case of tunneling between parallel surfaces, which would yield 
only 6% of the deviation. However, a considerable lowering and the strong spatial variation 
are in agreement with earlier measurements using a tungsten tip on a silicon (113) surface 
[5j]. A lowering was also observed in the first current- distance measurements at a voltage of 

n 

60 mV The variation of k at a subatomic scale can be explained with geometry effects 
in the orbital model for the tunneling process P, ^, Q| as shown above. At some adatoms 
of the silicon (lll)-(7 x 7) unit-cell, h shows a maximum as expected in the simple orbital 
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2 nm 




FIG. 4: Simultaneous measurement of topography (constant average tunneling current) (a), k 
(b), frequency shift A/ (c) and damping signal AE C per oscillation cycle (d). The image was 
acquired in the dynamic STM mode with a tungsten tip on the silicon (lll)-(7 x 7) surface and a 
qPlus sensor with a quality factor of Q = 2240. The surface unit cell is indicated with white 
diamonds. Sample bias: 2.0 V, average tunneling current I av = 0.8 nA, amplitude A = 0.1 nm. 

model. The lateral displacement of the maxima of R with respect to the topographic maxima 
is in very good agreement with the results of a model calculation for a s-tip on the silicon 
surface at the experimental conditions (8 ~ 20°). However, along the surface projection of 
A (vertical direction in the paper plane of Fig. @J, ^ decreases faster than in the model, 
probably caused by displacements of tip- and sample atoms at the small tip-sample distance. 
The variations in the maximal values of k between the four types of adatoms on the silicon 
(11 1)-(7 x 7) surface point to variations in the local work function. 

Partially, the lowering of the maximum value of k may be caused by the spatial extension 
of the tip atom or atom cluster. An extended tip results in a blurring of the image, compared 
to the theoretical images and should result in a lower value of k. In proximity to the sample 
the polynomial radial factor of the atomic basis functions can additionally cause a lowering 
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of k compared to k. Finally, the finite voltage causes changes in the barrier shape between 
tip and sample, which might also lower the value of k. 

The observed features may be a consequence of the structure of the surface wavefunctions 



which are probed in dynamic STM 



12j |. However, displacements of tip and sample 



atoms and the formation of chemical bonds also have to be considered because of the very 
small tip-sample-distance [3]. Figure 4 (c) shows that the frequency shift is less negative on 
top of the Si adatoms than between the adatoms. Because adatoms exert strong attractive 
forces on the tip before contact jl6l |. we conclude that the short-range forces acting between 
tip and sample in Fig. 4 are already repulsive. Bond formation can additionally cause 
variations in the decay constant images on an atomic scale. By which extent the measured 
variation is caused by mechanical deformations of the tip and/or sample, the formation of a 
chemical bond or by geometry effects solely because of the structure of the surface electron 
states can not be determined at this stage. 

In conclusion, we have presented simultaneous measurements of the tunneling current, 
the work function and the tip-sample forces with atomic resolution. Future improvements of 
the technique could be achieved by utilizing force sensors with even higher spring constants, 
such that smaller oscillation amplitudes become possible. Orienting the sensor oscillation 
exactly perpendicular to the surface would facilitate the interpretation of the data. We 
suggest that differences in the local barrier heights might help to identify the charge state 
or species of single atoms. The topographic corrugation can be increased dramatically using 
higher angular momentum tip states For atomic scale measurements with such 

tips, an enhanced corrugation of the work function images and higher harmonics of ln(///o) 
is expected. 
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